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ABSTRACT 

Using a systematic broad-band search technique, we have carried out a survey for large Lya nebulae 
(or Lya "blobs") at 2 < z < 3 within 8.5 square degrees of the NO AO Deep Wide-Field Survey 
(NDWFS) Bootes field, corresponding to a total survey comoving volume of ^10 8 h^ 3 Mpc 3 . Here, 
we present our spectroscopic observations of candidate giant Lya nebulae. Of 26 candidates targeted, 
5 were confirmed to have Lya emission at 1.7 < z < 2.7, four of which were new discoveries. The 
confirmed Lya nebulae span a range of Lya equivalent widths, colors, sizes, and line ratios, and most 
show spatially-extended continuum emission. The remaining candidates did not reveal any strong 
emission lines, but instead exhibit featureless, diffuse blue continuum spectra. Their nature remains 
mysterious, but we speculate that some of these might be Lya nebulae lying within the redshift desert 
(i.e., 1.2 < z < 1.6). Our spectroscopic follow-up confirms the power of using deep broad-band imaging 
to search for the bright end of the Lya nebula population across enormous comoving volumes. 
Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift — surveys 



1. INTRODUCTION 

Giant radio-quiet Lya nebulae (also know as "Lya 
blobs") provide a window into th e physics of ongoing 
massive galaxy f ormation 
Ivison et all 119981: 



2004; 
2006; 



Matsuda et al 



Ma 



Stei dgl et al.l 
120041: 



e.g. , iFrancis et al.l Il996 
' l200T 



„ jDev et al.l 1 2005: Sa ito et al.l 
jNilsson et al.l 120061: ISmith fc Jarvisl 120071 : 
Prescott et all 120091 : lYang et al.l I2009D . Rare systems 
found primarily in overdense regions ([P rescott et al.l 



20081: IMa tsuda et al. 2004, fisM EoO; ISaito et alJ l2Q0a 
Yang et al] I2009L l2010f h Lya nebulae are extremely 
luminous {Lhya °° 10 44 erg s _1 ) and are frequently 
associated with young, star- forming galaxy popu l ations 
or obscured A GN (e .g.. iBasu-Zych fc Scharl 1 2004: 
IMatsuda et all 120041: iDev et al.l _ l2005l : iGeach et al.l 
120071 120091 : Prescot t et al. | l2012bh . Theoretical 
and observational investigations into the power 
source behind Lya nebulae have painted a varied 
picture, present i ng ar gu ments for AGN p o wering 
(IChapman eQ] BOOllBasu-Zvch fc Scharfl 12001 
IGeach et al.l 12007 , 120091) . starburst-driven winds 
(iTaniguchi fc Shiova 120001 : iTaniguchi etHI 120011 : 
Matsuda et al. 2004) , spatially-extended st ar formation 
(Matsuda et al. 2007; Prescott et al. 2012b), cold accre- 
tion (iNilsson et al. 2006; S mith et al. 2008; Goe rdt et afl 
20101: iFaucher-Giguere et al.l 120101 : iRosd ahl & Bl aizotl 
2012 ) , or some c ombination t hereof ( Dey et all 
2005; Prescott etHI 120091 l2012ri iWebb et al.l 12009: 



Colbert et al.ll201lD . Larger samples of giant radio-quiet 
Lya nebulae drawn from unbiased surveys are needed in 
order to accurately measure the space density of these 
sources, to determine the emission mechanisms primarily 
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responsible for the Lya nebula class as a whole, and to 
understand their relationship to the more well-studied 
Lya haloes f ound around m any high redshift radio 
galaxies (e.g., [M cCarthy 199 3|, and references therein; 
van Oiik et all fl9 96: Weiding eret al.l 120051: iMilev et al 



2006: iBarrio et al.l 120081 : ISmith et alJ l2QQ9t iZirm et al 



2009). 



Increasing the sample of such a rare class of objects re- 
quires surveying large comoving volumes, but the stan- 
dard approach relying on narrow-band imaging is lim- 
ited by observational expense. Our complementary ap- 
proach was to design a new search technique using deep 
broad-band imaging ([Prescottl 120091 ). We employed this 
method to select a sample of Lya nebula candidates from 
the Boote s Field of the NOAO Deep Wide-Field Survey 
(NDWFS; Uannuzi fc Devi 119991 ). A full presentation of 
our survey algorithm as well as the selection function of 
this approach and the implications for the space den- 
sity of Lya nebulae c an be found in companion papers 
(|Prescott et al.l l2012ai 2013, in preparation; hereafter, 
Papers I and III, respectively). 

In this work (Paper II), we focus on the spectroscopic 
follow-up of Lya nebula candidates drawn from deep 
archival broad-band imaging of the NDWFS Bootes field. 
In Section [2j we give a brief outline of the survey de- 
sign and candidate sample. In Section [3j we discuss 
our spectroscopic observations and reductions. Section 
2] presents the 5 Lya nebulae discovered within this 8.5 
square degree survey and discusses our primary contam- 
inant sources, potentially an interesting population in 
their own right. Section [5] discusses the implications of 
this sample, and we conclude in Section [6] 

We assume the standard ACDM cosmology (Qm = 0.3, 
Qa — 0.7, h = 0.7); 1" corresponds to a physical scales of 
8.3-7.8 kpc for redshift s of z = 1. 2 — 2.9. All magnitudes 
are in the AB system (Oke 1971. 

2. SEARCH DESIGN 

We have designed an innovative search for Lya neb- 
ulae in deep broad-band imaging and applied it to the 
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NDWFS Bootes field (Paper I). Our search is most sen- 
sitive to the largest and brightest Lya nebulae because 
it leverages the deep blue (Bw) imaging of NDWFS to 
look for sources where bright Lya emission boosts the 
broad-band flux relative to the very dark sky. Thanks 
to the wide area («9 deg 2 ) of the NDWFS and the 
large width of the Bw filter (« 1275A, correspond- 
ing to Az ~ 1), our survey is able to probe an enor- 
mous comoving volume (^10 8 hf 3 Mpc 3 ) with archival 
data and significantly reduce the required observational 
overhead. Our survey is therefore complementary to 
smaller volume surveys for Lya ne bulae that rely on sen- 
sitive narrow-band imaging (e.g ., Matsuda et al. 2004; 
Saito et a l. 2006; Smith fc Jarvisll2007[ lYang et al.ll2009l 
2010; Matsu da et all 1201 If ). 

The search algorithm and candidate sample are dis- 
cussed in detail in Paper I. In brief, we selected giant Lya 
nebula candidates using wavelet analysis of the compact- 
source-subtracted Bw images. We selected a set of first 
and second priority candidates (Figure [T]) based on their 
Bw — R color, as measured using large 30 pix=7.7 // di- 
ameter apertures, an d wavelet size, as determ ined us- 
ing SourceExtr actor (jBertin fc Arnoutsl [1996) on the 
wavelet- deconvolved images. The final candidate sam- 
ple consisted of 39 first priority and 40 second priority 
sources over a search area of 8.5 square degrees. Both 
first and second priority samples contained candidates 
with diffuse morphologies (morphological category dif- 
fuse) as well as those that appear to be tight groupings 
of compact sources (morphological category group), as 
discussed in Paper I. In addition, we flagged 6 sources 
from outside these selection regions that showed promis- 
ing morphologies (third priority). 

3. SPECTROSCOPIC FOLLOW-UP 

In this section we describe the spectroscopic follow-up 
of the Lya nebula candidate sample. 

3.1. Observations & Reductions 

We targeted a total of 26 Lya nebula candidates (15 
first priority, 5 second priority, and all 6 third priority 
candidates). Longslit spectroscopic observations were 
obtained using the MMT and the Blue Channel Spec- 
trograph during 2007 May, 2008 April and June (Ta- 
ble [1]). The observations used the 300 1/mm grating with 
1//0/1//5 wide slits resulting in a resolution FWHM of 
2.6/3.4A and a wavelength range of AA w 3100 - 8320A. 

During our 2007 run, we chose slit orientations based 
primarily on the morphology of the candidate (i.e., 
aligned with the major axis of the emission as estimated 
from the .B^z-band morphology), and as a secondary cri- 
terion attempted to intersect a nearby bright reference 
object if possible. In practice, however, we found that 
the faintness of our candidates and the short duration 
of our spectroscopic exposures necessitated a positional 
reference. As a result, slit orientations during the 2008 
runs were chosen to always include a positional reference 
object, and consequently do not always trace the major 
axis of the Bw morphology. During these later runs, we 
dithered the target along the slit by « 5" between expo- 
sures to minimize the effect of any bad pixels. The full 
list of targeted candidates and the results of the spectro- 
scopic observations are given in Table O 



We reduced the spectroscopic data using IRAF0 We 
subtracted the overscan and bias, and applied a flat-field 
correction using normalized observations of the internal 
quartz flat-field lamps. Twilight flats were used to deter- 
mine the illumination correction for the science frames. 
We removed cosmic rays from the 2D sky-subtracted data 
using ^ap0 One-dimensional spectra were generated us- 
ing t he task apall a nd optimal variance- weighted extrac- 
tion (Valdes 1992); the spectral trace was determined 
using bright unresolved sources on the slit. We deter- 
mined the wavelength solution to an rms precision of 
^0.08-0.18A using HeArNe and HgCd comparison lamps, 
and then corrected the data for any slight systematic off- 
set using the night sky lines as a reference. The final 
wavelength calibration is accurate to ±0.42 A. The rel- 
ative flux calibration was based on observations of the 
standard stars BD+33 2642, BD+26 2606, BD+28 4211, 
Feige 34, and Wolf 13460 For each night, we applied a 
grey shift to compensate for any variable grey (i.e., inde- 
pendent of wavelength) extinction that may have affected 
a given standard star observation relative to one taken 
under better conditions. The sensitivity functions de- 
rived from individual standard star exposures were con- 
sistent to within <0.1 mag. 

Due to the faintness of the candidates and the fact that 
we are searching for luminous Lya nebulae at z « 2 — 3, 
the aim of our follow-up spectroscopic program was to 
look for strong, high equivalent width line emission. A 
single strong line in the blue can be identified as Lya 
rather than an unresolved [Oll]AA3726,3729 doublet (the 
only other possibility at these wavelengths) due to the 
fact that a detection of [Oil] would be accompanied by 
stronger detections of [O ill] AA4959, 5007 and Ha as well. 
Candidates with strong Lya emission at z ~ 2 — 3 in 
the 5v^-band are easily detectable with the MMT/Blue 
Channel down to a 5<r limit of « 1 — 7 x 10 _1T in 30 min- 
utes (assuming a FWHM 5 S = 12 A Lya line measured 
within a I'Ibxb" aperture). This corresponds to limit- 
ing line luminosities of ~ 0.3 — 5 x 10 42 erg s _1 , which 
are below the typical luminosities of giant Lya nebu- 
lae. However, continuum-only sources are much fainter 
and require longer integration times to yield high signal- 
to-noise ratio spectra, more than was available during 
our spectroscopic campaign. To allow us to target the 
largest number of candidates, we therefore carried out 
a quick reduction of the data in real-time and contin- 
ued integrating on each of the 26 targeted candidates up 
until the point where we could either confirm the pres- 
ence of a line or confirm the presence of continuum with 
no strong lines. The deeper spectroscopy necessary for 
studying the spectral properties of our confirmed Lya 
nebulae in detail as well as for detecting absorption fea- 
tures in continuum-only sources is left to future observa- 
tions with larger telescopes. 

4. RESULTS 

Of the 15 first priority and 5 second priority candi- 
dates targeted for spectroscopic follow-up, 4 first priority 

4 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

5 http://iraf.noao.edu/iraf/ftp/iraf/extern/xdimsum020627 

6 KPNO IRS Standard Star Manual 
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sources and 1 second priority source had confirmed Lya 
emission: we easily recovered the prev iously- discovered 
large Lya nebula at z « 2.66 (LABd05: l5ev et al,ll2QQ5h 
and discovered new, spatially-extended Lya nebulae at 
z w 1.67, z w 1.88, z w 2.14, and 2 w 2.27. In addition, 
we also targeted 6 third priority candidates that showed 
promising diffuse morphologies upon visual inspection. 
However, no Lya or [Oil] line emission was confirmed in 
any of the third priority candidates. In this section, we 
describe each of the confirmed sources in turn and then 
discuss the primary contaminants to our survey. 

4.1. Confirmed Lya Sources 

Figures [2][6] show the postage stamps, two-dimensional 
spectra, and one-dimensional spectra of the Lya sources; 
the measured properties are listed in Table [3l The 
spectral extraction apertures were chosen to maximize 
the signal-to-noise ratio of Lya. Redshifts were de- 
termined from the centroid of a Gaussian fit to the 
observed Lya line. No correction was included for 
Lya absorption, a potential source of bias in our red- 
shift estimates. The Bw sizes, isophotal areas, and 
surface brightnesses were measured above the median 
la surface brightness limit of the entire NDWFS sur- 
vey from the original Bw images using SourceExtrac- 
tor (det ect_minarea=5, dete ct_thresh=28.9 mag 
arcsec -2 ; iBertin fc Arnoutsl ll996). The Lya sizes were 
measured from the 2D spectra using SourceExtractor 
above the la surface brightness limit at the location of 
Lya (detect_minarea=5, detect_thresh^ 1 x 10" 18 
erg s _1 cm -2 A -1 arcsec -2 ; see Table [3]). The Bw sizes 
along the slit can underestimate the Lya sizes measured 
from the 2D spectra; in the case of PRG3, our deepest 
spectrum, the Bw size underestimates the Lya size by a 
factor of >1.3. 

Estimates for the total Lya isophotal area and total 
Lya luminosity are also given in Table [3j The approxi- 
mate total Lya isophotal area was estimated by correct- 
ing Ab w by a factor of ^ 2 , where Ab w is the isophotal 
area of the source on the Bw image and v is the ratio 
of the Lya and Bw sizes measured along the slit. The 
approximate total Lya luminosity was derived by scaling 
the Lya luminosity within the spectroscopic aperture by 
the geometric correction factor f geo = Ab w x v 2 / (uj x d) , 
where uj is the slit width and d is the spatial extent of the 
spectral extraction aperture. We stress, however, that in 
using area corrections based on broad-band imaging we 
are relying on the assumptions that the Bw emission 
roughly traces the Lya emission in the source and that 
the relative factor between the Bw and Lya sizes is the 
same throughout the object as it is along the slit; either 
of these assumptions may be violated in practice. 

In what follows, we discuss each confirmed Lya source 
in detail. Of the confirmed Lya nebulae, four were origi- 
nally categorized as having a diffuse morphology while 
one (PRG4) was categorized as having a GROUP mor- 
phology (see Paper I). 

4.1.1. PRG1 

PRG1 is a remarkable Lya n ebula (Figure [2j). As dis- 
cussed in Pre scott et al.l ((2009), it is the first example of 
a Lya nebula with strong, spatially-extended Hell emis- 
sion and weak metal lines, suggestive of a hard ioniz- 
ing continuum and potentially low metallicity gas. The 



Bw imaging shows a diffuse nebula and several compact 
sources, the brightest of which is located at the north- 
west edge of the nebula. Despite the strong Lya emission 
and large size (>78 kpc), PRG1 was selected as a sec- 
ond priority candidate because, at z ~ 1.67, Lya is at 
the edge of the optical window and not contained within 
the Bw band, giving the source a relatively red Bw — R 
color. Thanks to its diffuse blue continuum (95%) and 
Hell emission (5%), however, this source was still se- 
lected by our survey. When first discovered, this source 
was the lowest redshift Lya nebula known and the only 
one that shows strong spatially-extended Hell emission; 
PRG1 is therefore an ideal target for detailed study of 
the physical conditions and kinematics within Lya neb- 
ulae. Analysis of th e metallicity and sour ce of ionization 
in PRG1 is given in lPrescott et al.l (|2009[ ). and more de- 
tailed analysis using deep Keck/LRIS spectroscopy is in 
progress (Prescott et al., in preparation). 

4.1.2. PRG2 

PRG2 is a large Lya nebula at z ~ 2.27 with a 
roughly diamond-shaped morphology in the Bw image 
(Figure |3J). The identification of the strong line in the 
spectrum as Lya is secure based on the fact that no 
other lines are well-detected in the discovery spectrum 
and corroborated by weak detections at the positions of 
CivAA1548,1550, HeiiA1640, and Cm]A1909. In the case 
of [Oil] at lower redshift, we would have easily detected 
[Om] and Ha instead. The Lya nebula spans almost 
100 kpc, and at the southwestern corner there is a very 
blue compact source that appears to be a Lya-emitting 
galaxy from the spectrum. The redshift of this source 
is ideal for follow-up NIR spectroscopy as the rest-frame 
optical emission lines ([Oil], [Om], H/3, and Ha) will be 
observable in the J, H, and K bands. Continuum emis- 
sion is observed from two compact knots located at either 
end and from spatially-extended emission at fainter levels 
in between them. 

4.1.3. PRG3 

PRG3 is a Lya nebula at z ~ 2.14 (Figure [4j). It has a 
rather clumpy horseshoe-shaped morphology in the Bw 
imaging and spans ~74 kpc. The single strong line is 
identified as Lya rather than [O il] based on the fact that 
we do not see corresponding detections of [Om] and Ha. 
The spectrum shows spatially-extended continuum, but 
no other strong emission lines. 

4.1.4. PRG4 

PRG4 appears to be a candidate that was selected due 
to a close grouping of compact blue sources (Figure EJ). 
Due to the very blue color, it was flagged as a high 
priority target. At these wavelengths (blueward of the 
restframe wavelength of [Oil]), Lya is the only possible 
strong line. In addition, no other strong lines are seen in 
the spectrum. Although the Bw size of the full group- 
ing is roughly 7", the observed Lya at z « 1.89 is only 
marginally extended along the direction of the spectro- 
scopic slit (3.9", 33 kpc). The source may be larger in 
Lya: there is additional diffuse emission outside the slit 
that is visible to the southwest in the Bw imaging, but 
without further spectroscopy, we cannot determine if it 
is associated with coincident Lya emission. 
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4.1.5. LABd05 

LABd05 is the source that was the insp iration for 
our b road-band Lya nebula search (Figure [6j Dev et ah! 
20051). One of th e largest Lya nebulae known (> 100 kpc; 
Dev et al. 2005), it is located at z ~ 2.656. Our shallow 
MMT spectrum was taken at a slightly different posi- 
tion than the existing deeper spectroscopy from Keck 
but shows a hint of Hell emission and an emi ssion line 
at 5081A, both seen previously in the system ([Dev et al.l 
2005). The emission line at 5081A is thought to be 
Lya from a background interloper galaxy at z « 3.2, 
the compact source that is visible in the ground-based 
imaging and located at the western edge of the slit for 
this observation. Detailed study of ground-based data 
as well as high resolution imaging from HST showed 
that there are numerous compact galaxies, including a 
spectroscopically-confirmed Lyman break galaxy, within 
the sy stem that are offset spatially from the Lya nebula 
itself ([Dev et al.ll2"005l : iPrescott et al.ll2012bl ). The HST 
imaging demonstrated that the nebula contains diffuse 
restframe UV continuum emission, that the Lya emis- 
sion itself is smooth with a relatively round and disk-like 
morphology, and that the Hell e mission is spatially ex - 
tended by « 0.6 -I" (« 5-8 kpc: IFrescott et al.ll2012bD . 

4.2. Survey Contaminants 

The dominant contaminants in both the first and sec- 
ond priority spectroscopic samples are sources with spa- 
tially resolved blue continuum emission but no visible 
emission lines. Despite the lack of strong line emis- 
sion in the Bw band, our morphological broad-band 
search selected these sources either due to sufficiently 
extended, blue continuum emission or due to a close pro- 
jected grouping of blue galaxies. A few examples of these 
continuum-only sources are shown in Figure [71 

Without deeper spectroscopy, we can only speculate 
as to the nature of these continuum-only sources. The 
largest cases within the candidate sample (sources 1+2 
and 3; Paper I) are so spatially extended (~15-86 // , 
which at z ~ 1.2 — 2.9 would imply physical sizes of 
« 130 — 710 kpc in the continuum) and irregular in mor- 
phology that they are almost certainly located within 
the Galaxy, perhaps low surface brightness Galactic re- 
flection nebulae. Since low redshift (z < 1.2) blue star- 
forming populations or low surface brightness galaxies 
(LSBs) would be expected to show [Oil], [Olll], or Ha 
emission lines in our spectra, some fraction of the remain- 
ing continuum-only contaminants may in fact be galaxies 
or Lya nebulae in the redshift desert (1.2 < z < 1.6), 
for which Lya is blueward of the atmospheric cut-off 
(^obs ^ 3100A) but for which [On] has been redshifted 
past the red end of our MMT/Blue Channel spectra 
(A 6s ^ 8320 A). One of the Lya nebulae confirmed by 
our survey (PRG1, at z ~ 1.67) is in fact below the red- 
shift where Lya is covered by the 5^-band. Instead, 
this source was selected by our survey primarily due to 
blue continuum emission, and it was only thanks to the 
excellent blue sensitivity of MMT/Blue Channel that we 
were still able to detect the Lya emission at ^3250A. 
The case of PRG1 lends credence to the hypothesis that 
at least a fraction of the continuum-only "contaminant" 
sources are in fact Lya nebulae atl.2<z<1.6. 

At the same time, however, our expectation from Pa- 



per I was that Lya nebulae in the redshift desert should 
make up roughly 25% of the candidate sample, under the 
optimistic assumption that the Lya nebula number den- 
sity does not evolve significantly with redshift. In prac- 
tice, we found continuum-only detections represented a 
much larger fraction (75%) of the target spectroscopic 
sample, suggesting that this explanation may not be the 
full story. While the presence of continuum emission 
in the spectra does confirm that these continuum-only 
sources are indeed real astrophysical objects and not ar- 
tifacts within the NDWFS imaging, deeper ground-based 
optical spectroscopy or UV spectroscopy from space will 
be required to confirm their origin on a case by case basis. 
At this stage, their nature remains mysterious. 

5. DISCUSSION 
5.1. A Successful Broad-band Lya Nebula Survey 

This work is the first demonstration of the feasibility 
of conducting systematic surveys for large Lya nebulae 
using deep broad-band imaging datasets. The primary 
advantage of our unusual survey approach is the enor- 
mous comoving volume that can be surveyed using deep 
archival datasets. In addition, since this search technique 
is best used in the blue where the sky is dark, the re- 
sulting Lya nebula sample is weighted to lower redshifts 
(z < 3) where we have the opportunity to undertake de- 
tailed studies of their properties. The obvious trade-off is 
that our approach is not as sensitive to Lya nebulae that 
are intrinsically faint, low surface-brightness, or compact 
in morphology, as discussed in Papers I and III. Our 
search, therefore, provides a measurement of the bright 
end of the Lya nebula luminosity function, nicely com- 
plementing standard narrow-band surveys that probe to 
fainter luminosities. 

The success rate for finding sources with Lya emission 
was w 27% for first priority and « 20% for second prior- 
ity candidates. Therefore, if we were able to target all the 
Lya nebula candidates in our sample, we would expect 
to find a total of ~ 18 Lya nebulae (~ 10 and ~ 8 from 
the first and second priority set, respectively). While one 
of the goals of our broad-band survey for Lya nebulae is 
to place constraints on the space density of these rare 
objects, a robust estimate of the space density requires 
a detailed analysis of the selection function and is be- 
yond the scope of the present paper. Here, we briefly 
discuss our detection rate in the context of traditional 
narrowband Lya nebula surveys. 

Based on the result s of the na r rowba nd survey car- 
ried out at z ~ 2.3 by I Yang et al.l ([20091 ), in Paper I we 
estimated the expected number of Lya nebulae in our 
survey volume to be ~60-400, assu ming a 100% detec - 
tion rate, the same detection limit as I Yang et al.l ((2009), 
and a constant volume density as a function of redshift. 
Instead, we have confirmed 5 Lya nebulae, and scaling 
these results to the unobserved candidates, we expect to 
find only 18. While this estimate is extremely crude, 
it does suggest that the space density of the detected 
Lya nebulae in o ur sample is lower than that of the 
Yan g et all j2009) samp le. Possible reason s for this dif- 
ference are that (a) the lYang et all ((2009) narrow-band 
survey is more sensitive to fainter, and therefore less lu- 
mi nous, Lya nebulae than our broad-band survey; (b) 
the lYang et all (|2009l ) survey does not exclude Lya neb- 
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ulae with bright central sources whereas our survey does 
due to the natu re of the morpholo gical search algorithm; 
and/or (c) the I Yang et al.l ([2009) survey is more sensi- 
tive to cosmic variance than our larger volume survey. 
We defer a more detailed discussion of the space density 
of Lya nebulae implied by our survey to Paper III. 

5.2. Dispersion within the Lya Nebula Class 

The power of a systematic survey is the opportunity 
it provides to find out what is common among a class 
of objects and also what the dispersion in properties 
is among members of that class. The four large cases 
in our sample (PRG1, PRG2, PRG3, LABd05) span 
nearly an order of magnitude in total Lya luminosity 
(50 — 170 x 10 42 erg s _1 ), show a range of Lya equivalent 
widths (~ 50 — 260 A), and are at least 70 — 100 kpc in 
diameter. Morphologically, the four large Lya nebulae 
all show clumps and knots of emission in the broad-band 
imaging. The brightest compact knot in PRG1 is very 
red while that in PRG2 is remarkably blue. In addition, 
all four show what appears to be diffuse continuum emis- 
sion in the ground-based spectroscopy. This could either 
be due to many unresolved clumps, or due to a continuum 
component that is truly spatially extended. Analysis of 
HST imaging of one system (LABd05) lends support to 
the latter hypothesis, revealing that most of the contin- 
uum in thi s one source is unreso lved even at high resolu- 
tion (0.1": lPrescott et al,ll2Q12bh - In three cases (PRG1, 
PRG2, LABd05) there is evidence for emission in other 
lines (e.g., Civ, Hen, or Cm]). 

Given that diffuse continuum emission will have a 
larger impact on the observed broad-band color than line 
emission, one might ask if our survey is biased towards 
finding lower equivalent width sources than narrow-band 
surveys. In fact, however, Figure [8] shows that our survey 
uncovered Lya nebulae with rest-frame equivalent widths 
comparable to those of luminous Lya nebulae found us- 
ing standard narrow-band surveys but over a much larger 
redshift range. 

6. CONCLUSIONS 

We have carried out an innovative and economical sys- 
tematic search for large Lya nebulae using archival deep, 
broad-band data. While our technique is only sensitive 
to the largest and brightest Lya nebulae, it is able to 
probe enormous comoving volumes (~ 10 8 h^ 3 Mpc 3 ) 
using existing deep broad-band datasets. The details of 
our search algorithm, the selection function, and implied 
space density are discussed in Papers I and III of this 
series. In this paper (Paper II), we presented details 
of our spectroscopic follow-up of Lya nebula candidates. 
Within our ~8.5 square degree survey area and a redshift 
range of z ~ 1.6 — 2.9, we confirmed 4 new Lya nebulae 
and recovered 1 previously known case. The brightest 4 
Lya nebulae have Lya luminosities of ~ 5 — 17 x 10 43 erg 
s _1 and sizes of >70 kpc. Our broad-band search found 
Lya nebulae with large Lya luminosities and equiva- 
lent widths comparable to those found with narrow-band 
surveys, but revealed a new common theme: at least 
some large Lya nebulae show diffuse, spatially-extended 
continuum emission. The primary contaminants in our 
survey are sources that show nothing but blue contin- 
uum in the optical range, some of which we suspect 



may be galaxies or Lya nebulae located in the redshift 
desert. Deep continuum spectroscopy and comparisons 
to GALEX photometry will be required to confirm this 
claim. This work uncovered the first example of a giant 
Lya nebula at z < 2 and has demonstrated the feasibility 
of using deep broad-band datasets to efficiently locate lu- 
minous Lya nebulae within enormous comoving volumes. 
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UT Date 


Instrumental Resolution 21 
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1.0-1.2 Clear, high winds 

1.0- 1.2 Clear, high winds 
1.3-1.7 Clear, high winds 

1.0 Mostly clear 

1.2-1.9 Clear, high winds 
1.0 Clear 
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Table 2 

Spectroscopic Targets 



Candidate Name 


Right Ascension 


Declination Priority 


UT Date 


Exposure 


Class a 


Notes 




(hours) 


(degrees) 




Time (s) 







( 2 ) 


MTM H 7T710 

NDW1 b 


J 143006. 


,9+353437 


14: 


:30: 


:06. 


,864 


35: 


34: 


36. 


73 


3 


2007 May 20-21 


2400 


Continuum 


Galactic? 


(3) 


NDWl b 


T "1 A o o A r* 

J 142846. 


,2+330819 


14: 


:28: 


:46. 


,228 


33: 


:08: 


19. 


42 


3 


2007 May 20-21 


10800 


Continuum 


Galactic? 


(10) 


NDWFb 


J143411. 


,0+331731 


14: 


:34: 


:10, 


,975 


33: 


17: 


31. 


26 


1 


2007 May 20 


1800 


Lya 


LABd05 


/l A \ 

(14) 


NDWl b 


J143512. 


o i on 1 r\ r\ 

,3+351109 


14: 


:35: 


:12. 


,336 


35: 


11: 


:08. 


62 


2 


2008 Jun 8-9 


7200 


Lya 


PRG1 


(18) 


NDW1 b 


J 143222 


,8+324943 


14: 


:32: 


:22. 


,768 


32: 


49: 


42, 


,67 


2 


2008 Jun 8 


1800 


Continuum 


_ 


(24) 


NDWFS 


J142614. 


.7+344434 


14: 


:26: 


:14 


,714 


34: 


44: 


34, 


.22 


1 


2008 Jun 8 


3600 


Continuum 




(26) 


NDW1 b 


J 142622 


,9+351422 


14: 


:26: 


:22. 


,905 


35: 


14: 


22. 


.02 


1 


2008 Apr 3 


3600 


Lya 


PRG2 


(29) 


NDWl b 


J 142526. 


,3+335112 


14: 


:25: 


:26. 


,332 


33: 


:51: 


:12. 


,16 


1 


2008 Jun 9 


3600 


Continuum 




(31) 


NDWl b 


J 142547. 


.1+334454 


14: 


:25: 


:47. 


,126 


33: 


:44: 


:54. 


.13 


1 


2008 Apr 30 


2400 


Continuum 




(33) 


NDW1 b 


J 142714. 


o i o A o I r" p" 

,8+343155 


14: 


:27: 


:14 


,791 


34: 


:31: 


:54, 


,55 


2 


2008 Jun 8 


1800 


Continuum 






IN LJ VV -T O 


T1 ZL^I 98 




14 


•^1 
,o± 


•98 
. zo. 


,245 




9fi 

, ZD. 


O 1 . 


Q1 


2 


9007 Mn\r 99 
ZUU 1 Ivlciy — — 




I m "i~ l niin rn 
wOIlllIlU UIll 




(40) 


NDWFS 


J142653. 


.2+343855 


14: 


:26: 


:53, 


,172 


34: 


:38: 


:55, 


,39 


1 


2008 Apr 30 


3600 


Lya 


PRG4 


(44) 


NDWFS 


J142927 


,8+345906 


14: 


:29: 


:27, 


,837 


34: 


:59: 


:06. 


,14 


3 


2007 May 21 


3600 


Continuum 




(52) 


NDWFS 


J143706. 


.6+335653 


14: 


:37: 


:06. 


.588 


33: 


:56: 


:52, 


.65 


2 


2007 May 20-21 


4800 


Continuum 




(58) 


NDWFS 


J142516. 


,6+324335 


14: 


:25: 


:16 


.629 


32: 


:43: 


:35. 


.47 


1 


2008 Jun 8 


3600 


Continuum 




(59) 


NDWFS 


J143412. 


,7+332939 


14: 


:34: 


:12 


.722 


33: 


:29: 


:39, 


.19 


1 


2008 May 20 


10800 


Lyo; 


PRG3 


(65) 


NDWFS 


J143207 


,2+343101 


14: 


:32: 


:07. 


.224 


34: 


:31: 


:01. 


,34 


3 


2007 May 22 


3600 


Continuum 




(66) 


NDWFS 


J142539, 


,9+344959 


14: 


:25: 


:39 


,859 


34: 


:49: 


:59, 


.19 


1 


2008 Jun 8 


3600 


Continuum 




(70) 


NDWFS 


J142753. 


.8+341204 


14: 


:27: 


:53 


.762 


34: 


:12: 


:04, 


,10 


1 


2007 May 20-21 


8400 


Continuum 




(71) 


NDWFS 


J142600. 


,8+350252 


14: 


:26: 


:00. 


.842 


35: 


:02: 


:52, 


,36 


3 


2008 Apr 3 


3600 


Continuum 




(72) 


NDWFS 


J142643. 


.9+340937 


14: 


:26: 


:43 


.850 


34: 


:09: 


:36. 


.82 


1 


2008 Jun 9 


3600 


Continuum 




(73) 


NDWFS 


J142722. 


,4+345225 


14: 


:27: 


:22, 


.408 


34: 


:52: 


:24, 


.74 


1 


2008 Apr 3 


3600 


Continuum 




(74) 


NDWFS 


J142620. 


,0+340427 


14: 


:26: 


:19, 


.982 


34: 


:04: 


:27. 


,01 


1 


2008 Apr 30 


2400 


Continuum 




(80) 


NDWFS 


J142548 


,3+322957 


14: 


:25: 


:48. 


.283 


32: 


:29: 


:56. 


,58 


1 


2008 Jun 9 


3600 


Continuum 




(82) 


NDWFS 


J142449. 


,8+324743 


14: 


:24: 


:49 


.761 


32: 


:47: 


:42, 


.61 


1 


2008 Apr 30 


3600 


Continuum 




(85) 


NDWFS 


J142533. 


.0+343912 


14: 


:25: 


:32 


.966 


34: 


:39: 


:11, 


.95 


3 


2008 Apr 3 


3600 


Continuum 





Candidate ID numbers in the first column are the same as in Paper I. 
a Spectroscopic targets were classified as either showing "Lya" or "Continuum" emission. 
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Table 3 

Lya Nebula Measurements 



PRG1 PRG2 PRG3 PRG4 LABd05 



Aperture (arcsec) 1.5 x 5.04 1.5 x 7.84 1.0 x 5.60 1.5 x 1.68 1.0 x 4.48 

\Lya obs (A) 3249.61 ± 0.39 3971.41 ± 0.13 3813.28 ± 0.90 3511.23 ± 0.67 4444.99 ± 0.31 

Redshift 1.6731 ±0.0003 2.2668 ±0.0001 2.1368 ±0.0007 1.8883 ±0.0005 2.6564 ±0.0003 

F Lya (10- 17 erg s- 1 cm- 2 ) 44.1 ± 4.0 49.2 ± 1.1 10.2 ± 1.2 10.3 ± 1.2 19.0 ± 0.9 

L Lya (10 42 erg s- 1 ) 8.2 ± 0.7 19.3 ± 0.4 3.5 ± 0.4 2.6 ± 0.3 10.9 ± 0.5 

Lya EW rest (A) 257.1 ± 29.1 127.3 ± 6.3 47.1 ± 6.4 88.7 ± 11.8 115.5 ± 9.5 

Lya FWHM o6s (A) 9.19 ± 0.60 8.52 ± 0.19 23.36 ± 7.90 6.51 ± 0.89 15.44 ± 0.70 

Lya cr v (km s -1 ) 361.2 ± 23.7 273.9 ± 6.1 782.1 ± 264.3 236.8 ± 32.2 443.3 ± 20.0 

*CIVA1550 (10- 17 erg s- 1 cm- 2 ) a 2.1 ± 1.1 1.8 ± 0.8 < 0.8 < 0.8 < 5.4 

*HeiiAi640 (10 -17 erg s- 1 cm- 2 ) a 5.8 ± 1.0 1.8 ± 0.9 < 0.9 0.7 ± 0.5 1.4 ± 1.3 

^bniAi909 (10 -17 erg s- 1 cm- 2 ) a 4.8 ± 0.9 3.0 ± 1.1 < 1.5 < 1.4 < 1.8 



Bw Diameter Along Slit b (arcsec) 


8.96 


10.12 


6.72 


6.75 


9.32 


Bw Isophotal Area b (arcsec 2 ) 


40.9 


73.2 


45.3 


38.7 


54.4 


B\y Surface Brightness 13 (mag arcsec -2 ) 


27.2 


27.0 


26.8 


27.0 


27.0 


Lya Diameter Along Slit c (arcsec) 


9.24 


12.04 


8.96 


3.92 


8.96 


Lya Diameter Along Slit c (kpc) 


78.3 


99.0 


74.4 


33.0 


71.3 


Approximate Lya Isophotal Area d (arcsec 2 ) 


43.4 


103.7 


80.4 


>5.9 f 


50.3 


Approximate Total L^ yOL e (10 42 erg s -1 ) 


47.2 ± 4.3 


170.2 ± 3.7 


49.6 ± 6.1 


>2.6 f 


122.8 ± 6.0 



a Line flux upper limits are 2a values. 

b B\y sizes, isophotal areas, and surface brightnesses measured from the NDWFS imaging using SourceExtractor with a detection threshold of 
28.9 mag arcsec -2 , the median la Bw surface brightness limit of NDWFS. 

c Lya sizes measured from the 2D spectra using SourceExtractor with detection thresholds of [2.5, 1.0, 0.8, 2.1, 1.5] xlO -18 erg s -1 cm -2 A -1 
arcsec -2 , the respective la line surface brightness limits at the position of Lya. 

d Approximate Lya isophotoal area computed by correcting the Bw isophotal area Ab w by a factor of v 2 , where v is the ratio of the Lya and 
Bw diameters measured along the slit. 

e Approximate total Lya luminosity computed by scaling the Lya luminosity measured within the spectroscopic aperture by a geometric correction 
factor of fgeo = Ab w X u 2 /(uj x d), where Ab w is the isophotal area of the source on the Bw image, v is the ratio of the Lya and Bw diameters 
measured along the slit, uj is the slit width, and d is the spatial extent of the spectral extraction aperture. 

f As the Bw emission is not an accurate tracer of the Lya emission in PRG4, approximate luminosity and area estimates have been replaced 
with lower limits derived from the spectroscopic data alone. 
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Figure 1. Size vs. Bw — R color of Lya nebula candidates selected using the broad-band morphological search pipeline. The wavelet 
size corresponds to the size of the source in the wavelet power map, but does not indicate the true nebular size of the object (Paper I). All 
candidates are shown as filled black circles. The first and second priority selection regions indicated (blue dashed lines) contain 39 and 40 
Lya nebula candidates, respectively. The spectroscopic targets that showed Lya are indicated using red stars, and those that did not show 
Lya: are marke d with light blue circles . The large filled black circle with an error bar represents the typical color of low surface brightness 
galaxies (LSBs; Haberzettl et al. 2007), and the histogram (plotted on a linear scale) represents the distribution of Byy — R colors for field 
galaxies in NDWFS, demonstrating that the colors of our final Lya nebula candidates are substantially bluer than typical LSBs and field 
galaxies. 
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Figure 2. Imaging and MMT spectroscopic observations of PRG1, a Lya nebula at z ~ 1.67. The top row shows the optical Bw, R, and 
/ imaging along with the slit used for follow-up spectroscopy. The central panel contains the two-dimensional spectrum versus observed 
wavelength, and the bottom panel presents the one-dimensional spectrum extracted from a 1.5x5. 04" aperture. The error spectrum (red 
line) and B\y bandpass (blue line) are shown, and the positions of common emission lines are indicated. The one-dimensional Lya emission 
line profile is shown in the inse t panel. Spectroscopic data presented are from the second night of observations only (UT 2008 June 9; for 
details see Prescott et al. 2009). 
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Figure 3. Imaging and MMT spectroscopic observations of PRG2, a Lya nebula at z ~ 2.27. The top row shows the optical Bw, R, and 
/ imaging along with the slit used for follow-up spectroscopy. The central panel contains the two-dimensional spectrum versus observed 
wavelength, and the bottom panel presents the one-dimensional spectrum extracted from a 1. 5x7.84" aperture. The error spectrum (red 
line) and Bw and R bandpasses (blue lines) are shown, and the positions of common emission lines are indicated. The one-dimensional 
Lya: emission line profile is shown in the small inset panel. 
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Figure 4. Imaging and MMT spectroscopic observations of PRG3, a Lya nebula at z ~ 2.14. The top row shows the optical Bw, R, and 
/ imaging along with the slit used for follow-up spectroscopy. The central panel contains the two-dimensional spectrum versus observed 
wavelength, and the bottom panel presents the one-dimensional spectrum extracted from a 1.0X5.6" aperture. The error spectrum (red 
line) and B\y and R bandpasses (blue lines) are shown, and the positions of common emission lines are indicated. The one-dimensional 
Lya: emission line profile is shown in the small inset panel. 
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Figure 5. Imaging and MMT spectroscopic observations of PRG4, a Lya nebula at z ~ 1.89. The top row shows the optical Bw, R, and 
/ imaging along with the slit used for follow-up spectroscopy. The central panel contains the two-dimensional spectrum versus observed 
wavelength, and the bottom panel presents the one-dimensional spectrum extracted from a 1.5x1. 68" aperture. The error spectrum (red 
line) and B\y and R bandpasses (blue lines) are shown, and the positions of common emission lines are indicated. The one-dimensional Lyo: 
emission line profile is shown in the small inset panel. The size of PRG4 is unknown; the source is very compact along the spectroscopic 
slit, but additional diffuse emission that may or may not be Lyo; is visible to the south in the B\v imaging. 
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Figure 6. Imaging and MMT spectroscopic observations of LABd05, a previously-known Lya nebula at z « 2.656 (Dev et al. 2005). 
The top row shows the optical B\y, R, and I imaging along with the slit used for follow-up spectroscopy. The central panel contains 
the two-dimensional spectrum versus observed wavelength and the bottom panel presents the one-dimensional spectrum extracted from a 
1.0x4. 48 /; aperture. The error spectrum (red line) and Byy and R bandpasses (blue lines) are shown, and the positions of common emission 
lines are indicated. The one-dimensional Lya: emissio n line profile is shown in the small inset panel. The emission line at 5081A from a 
background galaxy at z « 3.2 is labeled for reference (|Dev et al. 2005). 
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Figure 7. Examples of the continuum-only contaminant sources discussed in the text fSection l4.2[) . Two-dimensional MMT/Blue Channel 
Spectrograph spectra are shown versus observed wavelength, with the target source centered at 0" along the spatial dimension. 
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Figure 8. Rest-frame equivalent width (EW re st , top panel) and spectroscopic redshift {z svec , bottom panel) versus total Lyq luminosity 
(LLya) for Ly« nebulae reported in the literature (Steidel et al. 2000; Nilsson et al. 2006; Smith & Jarvis 2007; Saito et al. 2008; Yang et al. 
2009) and those found using our broad-band search (red stars). Our survey succeeded in discovering Lya nebulae with rest-frame equivalent 
widths comparable to those found by previous narrow-band surveys, but covers a large redshift range using just a single broad-band filter. 



